Some exact results for the velocity of cracks propagating in non-linear elastic models 



O 

o 

(N 

Ph. 
< 



T. M. Guozden and E. A. Jagla 

Centra Atomico Bariloche, Comision Nacional de Energia Atomica (8400) Bariloche, Argentina 

We analyze a piece-wise linear elastic model for the propagation of a crack in a stripe geometry 
under mode III conditions, in the absence of dissipation. The model is continuous in the propagation 
direction and discrete in the perpendicular direction. The velocity of the crack is a function of the 
value of the applied strain. We find analytically the value of the propagation velocity close to the 
Griffith threshold, and close to the strain of uniform breakdown. Contrary to the case of perfectly 
harmonic behavior up to the fracture point, in the piece- wise linear elastic model the crack velocity 
is lower than the sound velocity, reaching this limiting value at the strain of uniform breakdown. 
We complement the analytical results with numerical simulations and find excellent agreement. 
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I. INTRODUCTION 

The velocity of a crack propagating in a brittle mate- 
rial is known to be related to the sound velocity in the 
material. This general statement can be qualitatively 
justified by noticing that a crack is a sort of elastic dis- 
turbance, although of course of extreme non-linear na- 
ture. Thus it is not surprising that its velocity is related 
to the velocity of propagation of small amplitude elastic 
deformations. However, when we want to be more pre- 
cise about the relation between crack velocity and sound 
velocity, difficulties appear. In text book treatments of 
linear elastic fracture mechanics, it is suggested that the 
maximum crack velocity is given by the Raleigh velocity 
^fl-Q This limit is expected to be achieved at large 
driving forces (i.e., large applied strain), since for low 
driving forces the discrete (atomic) nature of the mate- 
rial may reduce the velocity drastically (this is called the 
lattice trapping effect). Experimentally, an increase of 
the velocity with the applied strain is observed in gen- 
eral, however the limiting Raleigh velocity is almost never 
achieved. Q Microscopical observation of crack paths in 
different kinds of samples reveal one source of this dis- 
crepancy: at velocities roughly close to f/j/S a straight 
crack path destabilizes, becoming wandering, and gen- 
erating side branches at larger velocities. Some people 
have claimed ^4] that if this effect is taken into account 
(i.e., the true microscopic crack path length is larger than 
the apparent macroscopic length) then the classical pre- 
diction is verified. But this cannot be claimed to be al- 
ways the case.p Even restricting to cracks propagating 
is a stationary fashion along a perfectly linear path, a 
careful analysis reveals that crack propagation velocity 
cannot be determined independently of the microscopic 
details close to the crack tip.jj Ql This means that a 
purely macroscopic analysis using continuous approxima- 
tions for the material leaves the velocity of the crack un- 
determined. This is the reason why detailed models of 
the breaking phenomena at the atomic scale are neces- 
sary in determining crack velocities. 

A class of fully consistent models on which crack ve- 
locities can be calculated (albeit numerically) are lattice 
spring models where the material is represented by a set 



of point masses joined by springs. |6J These springs can 
break when some threshold deformation is reached giv- 
ing rise to cracks in the form of connected sets of broken 
springs. 

It has been recently established[3, II] that the propa- 
gation velocity in this kind of model crucially depends 
on the presence of anharmonicities of the springs. These 
anharmonicities are also called hyperelastic effects. The 
most spectacular case is that to hyperelastic stiffen- 
ing (i.e., springs becoming stiffer at large deformation), 
that can produce supersonic crack propagation, some- 
thing that had been considered impossible in classical 
treatments of fracture. However, the case of hypere- 
lastic softening is by far the expected most common 
case, since most decohesion potentials typically interpo- 
late smoothly between the weakly deformed material and 
the broken material, in which the elastic constants are 
formally zero. In this case, and even in the absence of 
other effects such as crack velocity oscillation or crack 
branching, hyperelastic softening produces a noticeable 
reduction of the crack velocity. 

Even in the relatively simple class of lattice spring 
models, quantitative predictions of crack velocity is elu- 
sive, since, as already stated, breaking of the material 
is a form of non-linear behavior, and it is typically very 
difficult to find exact results for non-linear models. The 
situation is even worse in the presence of hyperelasticity, 
which is an additional source of non-linear behavior. 

In this paper we show that taking a continuous approx- 
imation in the propagation direction in a class of lattice 
spring models, some exact results can be obtained for 
the crack velocity even in the presence of hyperelastic 
softening. These results shed light on the effect of hyper- 
elasticity on crack propagation, and serve as a starting 
point for other (most likely numerical) studies in more 
realistic models. 

We have also implemented the model numerically and 
compared the simulated results with the analytical ones, 
finding excellent agreement. 




FIG. 1: A sketch of the model studied: a set of 2N contin- 
uous non-linear elastic chains (represented by the continuous 
lines) are coupled through perfectly harmonic interactions if 
the vertical distance between chains is lower than a threshold 
value «i,fc. If this value is exceeded (this may occur only be- 
tween chains ui and —u\) the two chains decouple defining 
the crack. In the figure, intact springs are shadowed. The 
crack advances to the right as the system evolves. Boundary 
conditions are fixed displacements imposed at the top and 
bottom of the figure. 



II. THE MODEL 



We consider a lattice spring model in a stripe geom- 
etry with fixed displacement mode III boundary condi- 
tions. A continuous description is implemented in the 
propagating direction (along the stripe, chosen to be the 
X direction), whereas a discrete model is considered in 
the perpendicular (y) direction. Thus the model con- 
sists of a set of a fixed number {2N) of continuous elastic 
chains as depicted in Fig. ^ Taking into account the 
symmetry of the system, we solve the equations only for 
the upper half of it, in which each chain is labeled by a 
discrete index j, 1 < j < N. We consider the (out of 
plane) mode III displacement of the chains, that is noted 
Uj{x) for chain j. Chain j is coupled to the two neigh- 
bor chains by harmonic springs. These springs can break 
when the coordinate difference between chains is larger 
than a breaking threshold that we note Ubk- We will al- 
ways assume that the crack propagates in the middle of 
the stripe, i.e., between chains ui and —ui. Our aim is to 
find the stable propagation velocity of this crack. Chains 
UN and —UN are coupled to the lateral sides of the sys- 
tem, on which fixed displacements are applied. The sides 
of the system can be formally introduced as chains un+i 
(and — uat+i), with un+i{x) = {N + 1/2)5. This defines 
5 as the nominal displacement between adjacent chains in 
the system. Hyperelasticity comes from the assumption 
that the spring constant of a chain changes from a low 
stretching value fco when \duj{x) / dx\ < Uni, to a value 
/co7 when \duj{x) / dx\ > Uni- Thus Uni can be appropri- 
ately called the non-linear threshold of the chains. In the 
present paper we consider only the case of hyperelastic 
softening, namely 7 < 1. In short, the model is defined 



by the equation: 
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and p being the density of each chain. 

We want to obtain the solution to this equation when 
the external strain S is in between two limiting values. 
The lowest possible value for propagation corresponds to 
the Griffith's threshold Sq, at which the elastic energy 
available in the system ahead of the crack equals that 
stored in the broken springs behind the crack. An easy 
calculation shows that Sq = Ubk/V2N + 1. On the other 
hand the maximum external strain that can be supported 
by the system is the one that would break the system 
even in the absence of any pre-existent crack. Clearly 
this stress for uniform breaking Sjj is given by 61/ — Ubk ■ 

A few remarks are in order. Our model is obviously 
anisotropic, as there are continuous chains along the x 
direction, whereas the system is discrete in the y di- 
rection. Another source of anisotropy lies in the fact 
that hyperelastic softening is introduced only inside the 
chains, but not in the inter-chain springs. We have pre- 
viously indicatedpl that in fact it is the hyperelasticity 
in the propagation direction that drives the non-trivial 
evolution of the system. There is no point in introducing 
hyperelasticity in the interchain springs, as this has no 
important effect in the dynamics and complicates greatly 
the analytical treatment. Note also that chains are not 
allowed to break, it is only the vertical inter-chain springs 
that break. This forces the crack to remain in the center 
of the stripe and avoids effects such as crack branching. 

The wave velocity inside each chain is given by V^ = 
y^ko/p. In the highly stretched case, the spring constant 
changes in a factor of 7, so w e can define the stretched 
wave velocity I/J as FJ = \/koj/p. We will solve the 
model under the assumption that there is a stable prop- 
agation of a crack in the middle of the stripe, with a 
velocity V. As we will see, this velocity -if not zero- will 
never be larger than Vw, nor lower than VJ. 



III. SCALING PROPERTIES OF THE 
SOLUTION 

Before presenting the analytical results we have de- 
rived, it is interesting to indicate some constraints on 
the solution that can be obtained using scaling argu- 
ments only. Let us suppose we have obtained the solution 
Uj(x,t) corresponding to a given set of parameters Uni, 



Ubk, 7, and some applied stress S, and that this solution 
has a velocity V. It is then immediate to verify that 
a rescaled solution au is also a solution of the problem 
for an applied strain aS, with the same velocity V if the 
parameters Uhk and Uni are rescaled to autk, and auni- 
This means that the velocity can be written as a function 
of the combinations Ubk/uni, and 6/uni. 

A less trivial scaling can be obtained by changing a 
solution of the form u{x,t), to a new form w{x,t) = 
u(Ax, Bt), and finding A and B and new coefficients of 
the model for w{x,t) to be a solution. The calculation is 
straightforward, and we present only the result, that can 
be stated as the fact that the velocity of the crack should 
be of the form 
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(3) 
where / is an undetermined function of the arguments 
and where we also made use of the previously obtained 
scaling. 

Already from this scaling relation a very important re- 
sult can be obtained. This corresponds to the case in 
which there is no hyperelastic softening, i.e., in which 
the spring constants within the chains remain always 
equal to fcp. This can be formally obtained by letting 
Uni go to infinity. Then the r.h.s. of Eq. |j2Jl becomes 
^1 — 'yf{N, 0, 0). Since obviously in this limit the crack 
velocity cannot depend on 7, the only possibility is that 
f{N, 0,0) = 0, which implies V = Vw This means that 
in the absence of hyperelastic effects the crack velocity is 
equal to the wave velocity for any S > da- 
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with the additional condition to be satisfied at the crack 
tip: u(0) = Mbfc/2. 

Non-trivial solutions of this non- linear equation of mo- 
tion can be obtained by matching the solution in the dif- 
ferent sectors. It can be shown in general that the crack 
tip (located at a; = 0) must also be the point that sepa- 
rates low and high stretching regions, i.e., \du/ dx\ < Uni 
for a; > 0, and \du/dx\ > Uni for x < 0. In fact, ac- 
cording to the differential equation, when the non-linear 
threshold is reached, there is a change of sign in the 
pre-factor of d'^u/dx'^, that passes from [1 — (T^/T4j)^] 
to [7— (V/Vw)'^]. But d'^u/dx^ cannot change sign, oth- 
erwise the first derivative du/dx would be an extreme 
at that point, and this is inconsistent since we assumed 
\du/dx\ < Uni to the right and \du/dx\ > u^i to the left. 
Then the change of sign of the pre-factor of d^u/dx'^ 
must be compensated by a change of sign on the right 
hand side of the equation, and this is only possible at the 
point where the system is breaking and the right hand 
side changes from 3u — 3(5/2 to u — 3<5/2. This justifies 
our statement that exactly at the crack tip, the value of 
\du/dx\ is Uni- 

For a; > the solution of the differential equation has 
the form 
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where we have already used the constraint u(0) = Ubkl'2., 
and from the requirement \du/dx\x=o = Uni we get the 
velocity as 



IV. EXACT SOLUTION FOR N = 1 

We present here the exact analytic solution of the pre- 
vious non- linear model in the case A^ = 1. This will be 
a reference result for the further discussion of the more 
interesting cases with TV > 1. 

Assuming a stationary propagation, we introduce the 
stationary solution u(x) = ui{x — Vt,0), where V is the 
propagation velocity to be determined self-consistently. 
We choose the new reference system in such a way that 
the crack tip is located at x = 0. Thus, we must search for 
solutions of the piece- wise defined equation (we eliminate 
the tilde for simplicity): 
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Ubk 



Unl 



(9) 



Note that the velocity is independent of the value of 7, 
and that it is consistent with Eq. Q, as it must be. We 
also see explicitly that V = Vw ii Uni — > 00 . To check that 
this is a consistent solution, we must verify that there is 
a reasonable form of u{x) for a; < 0. 

When 6 is large enough, the solution for a; < consists 
of a concatenation of similar pieces of the form 
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This kind of solution is sketched in Fig. [3 When two 
pieces of this form are matched together, the derivative 
of u has a jump. A momentum conservation condition 
must be satisfied at those points. In fact, the integral 
of the force on an infinitesimal piece of chain during the 
time in which it passes through the singular point, must 




FIG. 2: Exact solutions for N = 1, ui,k/uni — 2 and S/u„i — 
1.31 in a) and S/u„i = 1.23 in b). Jumps of ^- in a) satisfy 
a momentum conservation condition given by Eq. Illlll . In b), 
the jump in -^ at x = xq is described by Eq. 1131 . In this 



case, only the portion of chains between x 
is in the hyperelastic regime 
andui(O) = ^. 



xo and x — 



In both cases -p- 



= —Unl 



be equal to the change of momentum. It is obtained that 
this conservation requires in fact that the derivative of 
u a jump. The value of the derivative is equal on both 
sides, and its value is 



\du/dx\ = Unlil - ^)/{iV/V^f - 7). 



(11) 



Note in Fig. [2ta) how the sort of triangular oscillation 
extends to x — > — oo, i.e., the excess of elastic energy 
in the system remains in the form of kinetic and elastic 
energy far behind the crack tip. 

When S is reduced, the amplitude of the oscillation 
described by Eqs. H1U|I and Hll|l reduces also, and it 
vanishes at a particular value of S. For S values lower 
than this, the previous solution is not valid. It turns out 
that in this case the solution for a; < is singular in our 
continuous system. A way to "regularize" this singular 
solution is to consider the continuous chain as the limit 
of a discrete chain of point masses. When the number of 
masses per unit length goes to infinity, we can describe 
the solution that is obtained in the following way [see 
Fig. [2Ib)]: there is a region for a;o < a; < in which the 
solution has the previous form given in Eq. IjlOl) . For 
X < xq, the solution returns to the linear regime, and is 
composed by a smooth part and a singular part. The sin- 
gular part is an oscillation that has an amplitude which 



goes to zero in the continuum limit, but a frequency that 
diverges in this limit, in such a way that it can carry (in 
the form of kinetic energy) the excess of elastic energy in 
the system. The regular part has the form 
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and the matching conditions at a;o for this regular 
part correspond to have continuity of the function, i.e, 
m(xq") = u{xq), and a jump in the derivative given by 
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(13) 
which is obtained using the same kind of conservation 
arguments that led to Eq. (|II(I . In Section VI we will 
show results of numerical simulations that confirm and 
clarify further this behavior. 

We see in the solution for x < (Eq. I1U|) that when 
V -^ V2 the frequency of the oscillation for a; < di- 
verges. In fact, if according to Eq. |^the velocity would 
be lower than V.^, then neither the solutions given by 
Eq. H10|l nor Eq. H12|l exist. It can be shown that in this 
case the velocity of crack propagation is actually V — V2 ■ 
This regime is not particularly interesting to us, and from 
now on we will always assume to have chosen values of 7 
such that V >V;^. 

Note that according to Eq. jnj the crack velocity at 
the Griffith's threshold Sg — Ubk/VS^s finite iiutk/uni < 
2/(a/3 — 1), and is given by 
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This is in contrast to cases in which the system is discrete 
in the direction of crack propagation. In those cases, 
due to lattice trapping effects the crack cannot propagate 
if the Griffith's threshold is not overpassed by a finite 
amount. 

The values in Eq. ^ for the velocity as a function of 6 
must be compared with the result that is obtained in the 
absence of hyperelastic effect, namely V{6) = I^. The 
non-trivial result contained in Eq. JHJ is a consequence of 
hyperelasticity in the system. We will see that the same 
qualitative effects exist in the more interesting cases with 
N >1. 



V. EXACT RESULTS FOR A^ > 1 

The previous case A^ = 1 is a good starting point in 
which the analytical solution can be worked out in full 
detail. But obviously, if we are interested in modeling a 
macroscopic system we should study the case of a large 
number of chains. 

For A^ > 1 the exact value of the velocity for arbitrary 
d cannot be obtained, in general. However, we can pro- 
vide exact results in some neighborhood of Sg and Su. 



Consider first the case S ^ 5g- Let us concentrate on 
one half (the upper one) of the system, since the other 
is symmetric. Sufficiently close to the Griffith's thresh- 
old, only one chain (the one adjacent to the crack) will 
enter the hyperelastic region. In this regime the prob- 
lem can be separated in three sectors as shown in Fig|3 
We match the solutions, requiring continuity of the func- 
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FIG. 3: A case with N — 5 and 5 close to Sq- For clarity, we 
only plot the upper half of chains, as the others are symmetric. 
Only the chain adjacent to the crack, and only in region II, 
explores the nonlinear regime. As (5 — » (5g, 2;o — > 0. 

tion and derivative of Uj (x) , except for the derivative of 
ui{x), in which (as in the A^ = 1 case) a discontinuity 
of the derivative of the form H13|l exists between regions 
II and III. The solution obtained will be valid as long as 
no chain other than the first enters the non-linear regime 
and ui{x) < U2{x). The width of zone II is determined as 
part of the solution. The problem stands as a system of 
AN + 1 nonlinear algebraic equations, which we solve to 
any desired accuracy through an iterative method. The 
results for the velocity are plotted in Fig. 0]as a function 
oiS/Sc, and in Fig. [31as a function oiS/Su- We plot data 
in the full range in which the method is reliable and gives 
the exact value of the velocity. As it can be observed, the 
velocity is only weakly dependent on 7 (always assuming 
V^<V). 

It is interesting to observe from Fig 01 that even in 
the limit N —f 00 our method provides the solution in 
a finite range of S. This means that in all this interval, 
for a system of infinite chains, there is a single one that 
explores the hyperelastic regime, and is responsible for 
the full reduction of the velocity from Vw to the actual 
value. 

In the limit (5 — > (5gj we have xq —^ 0, and region II 
shrinks to zero. In this limit we obtain the exact values 
of the velocity and its derivative with respect to S by 
solving a linear system of 2N algebraic equations. Both 
V{Sg) and ^L turn out to be independent of 7. From 
this independence and the scaling form Eq. ||2Jl, we can 
conclude that 
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FIG. 4: Analytical results for the velocity vs S/5g, for 
Ubk/uni = 4/3. It is seen that the dependence on 7 is very 
weak. Through all the range in which we plot the data only 
chain ui (and — ui) explores the hyperelastic regime. 
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FIG. 5; Same as in the previous figure, plotted as a function 
oiS/Su- 



and 



Vi - (v/v^y 



dS 



Sg,N 



F2iN) 

Unl 



(16) 



The values of the functions Fi and F2 for different N 
are shown in Fig. |H1 The limiting value Fi{N -^ 00) = 
1/2 can be obtained analytically through an appropriate 
analysis of the equations in this limit. Extrapolation of 
the finite N exact values of F2{N) suggests also that 

limAr_^oo [F2{N)/ y^N/2] = 1, but we have not verified 
it analytically. 

We can then write: 
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We emphasize again the reduction of the velocity from 
the value Vw due to the hyperelastic effect. This effect 
disappears if Uni -^ 00. 

A second limiting case can be solved analytically, and 
that is the asymptotic form of the velocity very close to 
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the velocity can be written as 
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FIG. 6: Plot of the functions Fi(iV) and Fi^N) [see Eq2. 
itTCl and ifra i. 



Sij. The analysis is based again in matching the solu- 
tions of the piece-wise hnear equation. The situation is 
sketched in Fig[7| In sector I (a; > 0) all chains are in 
the linear regime. In sectors II, III, etc., chains enter 
the nonlinear regime sequentially, starting from the one 
adjacent to the crack. 




xo 

FIG. 7: Sketch of the configuration for S close to 5u ■ In sectors 
I, II, 111, etc., we have 0, 1, 2, etc. chains in the hyperelastic 
regime (dotted lines). The velocity for S close to 5u can be 
obtained analyzing sectors 1 and II only, as explained in the 
text. 



When S -^ 6u (and V/Vu, — > 1), it can be seen that 
the solution of the equations for sector II are N — 1 ex- 
ponential modes with a diverging decaying constant, and 
a trigonometric mode with finite frequency. Taking into 
account that the width of region II (namely |a;o|) remains 
finite even for V -^ Vw, 'we can neglect exponential modes 
that grow toward negative x. This allows us to obtain 
the velocity in this limit by solving a system of 2N linear 
equations, matching the solutions between regions I and 
II only. 

The result we obtain is that to lowest order in 6 — Su, 




S-Si 

Unl 



■FsiN) (19) 



S^St, 



where F^ is another A'^-dependent dimensionless func- 
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FIG. 8: Plot of F3{N) (see Eq. [THJ. The data follow a loga- 
rithmic trend, vanishing for A'^ — *■ oo as 1/lnA'^ as the fitting 
(dotted line) shows. 

tion. Note that this result is again independent of 7, and 
consistent with the general expression in Eq. ||2Jl. Values 
of F3{N) are plotted in Fig|Sl By analyzing in more de- 
tail the N ^' 00 limit, it can be shown that F3 goes to 
zero as l/ln(A^). 



VI. COMPARISON WITH NUMERICAL 
RESULTS 

Although the main results of our work are the analyti- 
cal findings of the previous section, we include here some 
results of numerical simulations for two reasons: first of 
all, some of the results of the previous section are not 
fully intuitive, and then we think it is clarifying to check 
them against a numerical simulation. Secondly, numer- 
ical simulations can be done in the full range between 
d = do and S = du^ filling the gap between the two ana- 
lytical limits. 

In our numerical simulations we are forced to consider 
a system that is discrete also in the x direction. We 
do this by introducing softer springs in the x direction 
than in the y direction. The ratio between horizontal 
and vertical spring constants will be noted k, and it is 
a measure of the degree of anisotropy of the lattice. For 
K = 1 the lattice is isotropic, whereas for k ^ we 
recover the continuous limit of the analytical treatments. 
As we will see, keeping a finite but small k is also an 
appropriate form of regularizing the singular results that 
may appear in the continuum case. 

In Figs. 1^ and [TUl we present results for A^ = 1. They 
compare very well with the analytical results of the pre- 
vious Section. Note in particular in Fig. Elb) the oscil- 
lation of du/dx for x < xq. This represents the singular 
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FIG. 9: a): Analytical (line) and numerical results (points) 
for N — 1, S/uni = 1.23 and Ubk/uni = 2 (k = 1/1600 in 
the numerics). For clarity, only one every four points of the 
simulated system is shown, b) The plot of du/dx, from the 
simulations. We calculate the derivative in the discrete system 



for chain j as du/dx ; 



/k [Ui+l^j — Ui^. 



, where the subindex i 



is the discretization along the x axis. Note that the only piece 
of chain in the non- linear regime is that with xq < a; < 0. The 
oscillation for x < xo carries (in the form of kinetic energy) 
the excess of elastic energy that is present in the system. 



behavior of the analytical solution that we have discussed 
previously. This oscillation is seen in the plot of du/dx, 
but is hardly visible in the plot of u{x) itself, as its am- 
plitude goes to zero with k. In the case of Fig. ^]we see 
how the abrupt change on the derivative of the analytical 
solution is very well reproduced in the numerical simula- 
tions, supporting the prescription given by Eq. Hll|) . In 
Fig. II II we show superimposed analytical and numerical 
results for N = 20, for a value of 6 in which a single chain 
explores the non-linear regime. Again the agreement is 
very good. Note also in this case the oscillation in du/dx 
for the first chain behind the crack. 

The numerical results for crack velocities are plotted 
on top of the analytical results in Fig. E| We see that 
they agree very well with the analytical values, providing 
a link between the 6 ^ 5g and 5 ^ 6u regions. 



VII. THE CONTINUOUS LIMIT: A^ ^ oo 

From the finite N results of the previous sections we 
can try to obtain the behavior of a 'macroscopic' material 
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FIG. 10: a): Analytical (line) and numerical results (points) 
for N = 1 with 5/uni = 1.53 and Ubk/uni = 2 (« = 1/1600 in 
the numerics). For clarity, only one every five points of the 
simulated system is shown. The solution behaves as described 
in Eq. dlj. The numerical results differ at the left border 
of the system, because the system is finite in the numerical 
simulation. 
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FIG. 11: a): Analytical (lines) and numerical results (points) 
for N — 20, 5/uni — 0.33 and utk/uni = 2 (k = 1/625 in 
the numerics). For clarity, only one every four points of the 
simulated system is shown. At this strain value only the first 
chain enters the non-linear regime, b) Numerical values of ^, 
where the singular oscillation behind the crack is observable. 



by studying the limit A^ — > oo. We must keep in mind 
however that the results we are about to discuss still de- 
pend strongly on the microscopic details of the model, 
other microscopic realization giving rise probably to dif- 
ferent macroscopic behavior. In fact, we have already em- 
phasized that the fracture of a macroscopic object cannot 
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FIG. 12: Numerical results for different number of stripes 
with 7 = and Ubk/uni = 2. Analytical results close to 5g 
are plotted as dashed lines, whereas the analytical asymptotic 
slopes at (5(7 [Ea. ll9ll 1 are shown by continuous lines. The 
agreement is seen to be excellent, considering the discreteness 
of the numerical model. Note that for the parameters chosen 
the velocity at 5g goes to zero when N ^> oo [Eq. II17II 1. 




FIG. 13: Same numerical results as in Fig. 1121 for different 
number of stripes with 7 = and Ubk/uni — 2, but plotted 
as a function of S/Sq- The numerical results suggest that, 
given S/Sa, the velocity is a decreasing function of A'^. This 
evidence allows us to plot an upper bound (continuous line) 
for the results in the N ^ oo case (see the text for details). 



leading analytical form can be obtained as 



be described completely in terms of a continuum descrip- 
tion, since microscopic details of the breaking process at 
the crack tip are always relevant. In any case, taking the 
N —f oo limit in our model provides us with one possi- 
ble realization of a continuum system that we want to 
analyze now. 

The most important quantity we can analyze in the 
large N limit is the dependence of the crack velocity on 
the normalized strain S/Sq. This is in fact a directly ac- 
cessible experimental quantity. We already have at hand 
some analytical results about this (see Eqs. 1171 and [TH|l . 
We know the value of V at Sq (which is strictly lower 
than Vw), and that of dV/d{S/SG) at S — Sq (which is fi- 
nite). We also know that eventually V reaches the value 
Viu for sufficiently large 5/6o- But the extremely slow 
decay with N of dV/dS at d ^ Su (Eq. lO) allows 
to infer that V{S/Sc) will reach the value Vw also very 
slowly. We present here a non-rigorous argument, which 
we think reproduces the right tendency. First of all note 
from the numerical results of Fig. E|that the asymptotic 
value of V close to Sjj is a reasonable upper bound for 
the velocity, i.e. [see Eq. ((T^ ]. 
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On the other hand, when plotted as a function of 5/5g 
(Fig. I13() . the numerical results suggest that, for a fixed 
value of 5/6g, the velocity is a decreasing function of iV. 
We think this is rigorous result, although we do not have 
a proof of it. Accepting this statement as valid, we can 
obtain an upper bound for the velocity in the N ^^ oo 
limit by maximizing the right hand side of Eq. (|20(l with 
respect to N for each value of S/6g- The result is plotted 
in Fig. ^1 as a continuous line. For large S/6g, the 
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(21) 
where a is a numerical constant. This is in fact an ex- 
tremely slow convergence to the limiting value Vw We 
think this is a very important result. It tells that strictly 
speaking, hyperelastic softening does not reduce the lim- 
iting value of the velocity for sufficiently large applied 
strain. However, in view of its extremely slow conver- 
gence to this limit, from a practical point of view we can 
say that hyperelastic softening produces an appreciable 
reduction of the limiting velocity. In particular, we see 
that this reduction of the velocity does not exist when 
hyperelasticity is absent (namely, for Uni -^ oo). 



VIII. DISCUSSION AND CONCLUSIONS 

We have analyzed the effect of hyperelastic softening in 
a model of crack propagation in a stripe geometry under 
mode III fixed displacement boundary conditions. The 
model is continuous in the propagation direction and has 
a finite number of chains in the perpendicular direction. 
The two central chains of the stripe can decouple when 
they separate more than a critical distance Ubk , generat- 
ing a crack in the model. 

In the case in which the chains are perfectly harmonic 
the velocity of crack propagation is equal to the wave ve- 
locity Vw in the full range of external strain 5 between the 
Griffith's threshold 5g and the strain of uniform break- 
down (5(7. 

We have studied how this result is affected by the in- 
clusion of hyperelastic softening in the chains, namely, a 
softening of the spring constant of the chains when the 
stretching is greater than a threshold value. We have pro- 
vided analytical results in some cases, and complemented 



them with numerical simulations. 

For the case of a single chain the full analytical solu- 
tion has been worked out. It is clearly seen already in 
this simple case that hyperelastic softening reduces the 
velocity from the harmonic case. Now the velocity has 
a non-trivial dependence on 6, and becomes equal to Vw 
only at Su- 

We have given the analytical solution for the veloc- 
ity in neighborhoods oi 5 — Sg and S — Sjj- The main 
results in this case are the following. The crack veloc- 
ity at 5g is strictly lower than V^. It decreases as a 
function of the number of chains but may well be finite 
in the N ^ oo limit for some range of the parameters 
of the model. There is a finite range of S/dc in which 
only the chain adjacent to the crack enters the hypere- 
lastic regime. This range remains finite for large N. This 



means that our analytical treatment provides the exact 
value of the velocity in a finite range around S/Sg = 1 
even in the N ^ oo case. 

The crack velocity tends always to Vw when S ^ Sjj. 
In the large N limit this can be stated as the fact that V 
tends to Vw for S/Sg -^ oo. However, ours estimations 
show that this convergence is very slow, namely like ^ 
l/ln\S/6G). 
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